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t  >  5  ms,  a  second  peak  emerges  within  our  de 
tectable  range  near  1/pm  <k<  1 .5/pm.  The  signal 
to  noise  ratios  for  the  peaks  are  statistically  limited 
to  1  to  2  for  images  taken  at  a  fixed  hold  time. 
Combining  data  points  at  all  times,  one  recovers 
clear  oscillations  as  observed  in  Fig.  2B.  The 
multiple  peaks  and  troughs  in  momentum  space 
result  from  the  interference  of  sound  waves  with 
different  wavelengths  and  are  the  key  features  of 
Sakharov  oscillations  in  the  angular  spectrum  of 
the  CMB  radiation.  Moreover,  moving  acoustic 
peaks  as  time  evolves  suggests  that  the  cor 
relations  are  spreading  out  at  a  finite  speed  of 
sound.  For  the  quench  up  experiments,  the  os 
cillations  are  less  clear,  owing  to  the  limited 
imaging  resolution  and  the  sample  size.  None 
theless,  we  clearly  observe  that  correlations  are 
suppressed  at  smaller  k  as  time  evolves.  Spread 
ing  of  correlations  along  the  acoustic  “light  cone” 
was  also  reported  in  a  quenched  atomic  Mott 
insulator  (26). 

The  oscillations  in  the  time  and  space  domains 
can  be  understood  as  the  interference  of  counter 
propagating  phonon  pairs  b\ke±lkr  ie(k)x/* *9  where 
is  the  creation  operator  of  phonons  with 
momentum  ±k  and  energy  e(k).  Assuming  that  pho 
nons  are  long  lived,  the  structure  factor  follows  a 
simple  form  (10,  20),  S(k, t)  S(k, 0)  +  2 Bk  cos(|) 
sin2(7i/x),  where  Bk  >  0  is  the  oscillation  amplitude, 
the  oscillation  frequency  is  given  by  file)  2e(k)/h, 
and  the  relative  phase  is  4>  0  for  quench  up 

experiments  and  4>  n  for  quench  downs  (20). 
Adopting  sinusoidal  fits  to  the  measured  data  in 
Fig.  2,  we  determine  the  oscillation  frequencies 
for  various  k  and  gf.  The  result  in  turn  provides 
good  fits  to  the  data  shown  in  Fig.  3  (20). 

We  summarize  five  sets  of  measurements 
on  quenched  superfluids  and  compare  the  os 
cillation  frequencies  to  the  predictions  from  the 

Bogoliubov  theory:  e  (k)  =  vhk  yj  \  +  k2  /  4,  where 
v  =  h  yjngf  /m  is  the  speed  of  sound  and  £,  h/mv 
is  the  healing  length.  Plotting  the  frequencies  in 


the  dimensionless  units  based  on  calculated  v  and 

the  data  in  all  three  quench  down  measurements 
collapse  to  a  single  curve,  which  agrees  with  the 
Bogoliubov  theory  (Fig.  4B).  Quench  up  experi 
ments,  however,  show  consistently  65%  higher 
frequencies  than  the  prediction  (Fig.  4A).  The 
discrepancy  likely  results  from  the  effects  of  strong 
initial  density  fluctuations  in  quench  up  experi 
ments  and  finite  temperature  (20). 

The  space  time  dependence  of  the  density  cor 
relations  in  quenched  superfluids  not  only  supports 
the  picture  of  Sakharov  acoustic  oscillations,  but 
also  suggests  that  Sakharov  oscillation  is  a  uni 
versal  feature  of  a  quenched  hydrodynamic  sys 
tern  far  from  thermal  equilibrium.  The  time  and 
spatial  scales  of  the  oscillations  are  related  by  the 
speed  of  sound  v,  and  the  observed  propagation 
of  Sakharov  acoustic  peaks  is  consistent  with  the 
sonic  horizon  k  u/vt  within  which  the  waves 
can  interfere.  In  particular,  good  agreement  with 
the  Bogoliubov  theory  is  found  when  the  inter 
action  is  quenched  to  smaller  values,  and  the 
growing  multipeak  structure  in  the  density  noise 
power  spectrum  resembles  that  in  the  CMB  ra 
diation  spectrum.  The  observation  of  Sakharov 
oscillations  in  a  quenched  atomic  superfluid  raises 
prospects  of  investigating  universal  hydrodynamic 
phenomena  in  a  laboratory  setting,  as  well  as  of 
extending  the  current  experimental  tools  to  address 
exciting  topics  in  cosmology  and  gravitational 
physics  such  as  Hawking  radiation  (13)  or  Unruh 
effect  (27). 
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and  tests  of  fundamental  physics  (1,  2).  Under 
development  worldwide,  optical  clocks  oscillate 
105  times  faster  than  their  microwave  predeces 
sors,  dividing  time  into  finer  intervals  (3).  Although 
microwave  clocks  such  as  the  Cs  fountain  have 
demonstrated  time  and  frequency  measure 
ments  of  a  few  parts  in  10 16  (4,  5),  optical  clocks 
now  measure  with  a  precision  of  1  part  in  1017 
(6  9). 

A  clock’s  instability  specifies  how  its  tick 
ing  fluctuates  over  time,  a  characteristic  generally 
quantified  by  the  Allan  deviation  (10).  No  time 
or  frequency  standard  can  make  measurements 
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Fig.  1.  Experimental  realization  of  the  Yb  optical  lattice  clocks.  (A)  Laser 
light  at  578  nm  is  prestabiized  to  an  isolated,  high-finesse  optical  cavity  using 
Pound-Orever-Hall  detection  and  electronic  feedback  to  an  AGM  and  laser 
piezoelectric  transducer.  Fibers  deliver  stabilized  laser  light  to  the  Yb-1  and  Yb-2 
systems.  Resonance  with  the  atomic  transition  is  detected  by  observing  atomic 
fluorescence  collected  onto  a  PMT.  The  fluorescence  signal  is  digitized  and  pro¬ 
cessed  by  a  microcontroller  unit  (MCU),  which  computes  a  correction  frequency, 


/12( f).  This  correction  frequency  is  applied  to  the  relevant  AGM  by  way  of  a  dfrect 
digital  synthesizer  (DOS)  and  locks  the  laser  frequency  onto  resonance  with  the 
clock  transition.  (S)  Relevant  Yb  atomic  energy  levels  and  transition^  including 
laser-cooling  transitions  (399  and  556  nm),  the  dock  transition  (578  nm),  and  the 
optical-pumping  transition  used  for  excited-state  detection  (1388  nm).  (0  A  single¬ 
scan,  normalized  excitation  spectrum  of  the  dock  transition  in  171Yb  with 
140-ms  Rabi  spectroscopy  tine;  the  red  line  is  a  free-parameter  an c2  fuiction  ft 


Fig.  2.  Frequency  com¬ 
parison  between  the  Yb 
optical  lattice  docks.  (A) 
Correction  frequencies, 
are  shown  in  red 
and  blade  Dominant  LO 
fluctiBtions  are  due  to  the 
tavrty  and  are,  thus,  com¬ 
mon  to  the  atomic  systems. 
(B)  Frequency  difference, 
/2(f)  - /i(t),  between  the 
two  Yb  clock  systems  for 
a  5COO-S  friterval.  (C)  Data 
set  /2(ti  -  /i(f)  over  a 
90,000-s  inteval  Gaps  rep¬ 
resent  data  rejected  before 
data  analysis  due  to  servo 
unlocks.  (D)  Histogram  of 
all  data  and  a  Gaussian  ft 
{Xf  =  0.9996). 
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better  than  the  statistical  precision  set  by  its  in 
stability.  Further,  a  dock's  systematic  uncertainty 
is  often  constrained  by  its  long  term  instability. 
For  these  reasons,  and  because  many  timing  ap 
plications  require  only  exquisite  instability,  the 
instability  represents  perhaps  the  most  important 
property  of  an  atomic  standard. 


In  pursuit  of  lower  instability,  an  optical  lat 
bee  clock  uses  a  stabilized  laser  referenced  to 
many  (103  to  106)  alkali  ne  earth  (o  r  similar)  atoms 
confined  in  an  optical  standing  wave.  Alignment 
of  the  clock  interrogation  laser  along  the  direc¬ 
tion  of  tight  lattice  confinement  eliminates  most 
Doppler  and  motional  effects  while  probing  the 
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Fig.  3.  Measured  instability  of  a  Yb  optical  lattice  clock.  Total  Allan  deviation  of  a  single  Yb 
clock,  - /tWWs  (red  rirdes),  and  its  white-frequency-nose  asymptote  of  3.2  x  MT16/\/t  (red 
solid  Hie).  The  blue  dashed  line  represents  the  etfimated  combined  instability  contributions  from  the  Dick 
effect  (1,4  x  10 ~uly/z)  and  QPN  (1  x  1(T16/vM;  the  shaded  region  denotes  uncertainty  h  thee 
estimates  Error  bars  indicate  la  confidence  intervals 


Averaging  time  t(s) 


Fig.  4.  Calculated  instability  limits  toward  the  goal  of  1  x  10_1S  in  100  s.  The  Dick  limit  (red 
dotted  line)  is  reduced  by  using  a  hypothetical  LO  four  times  as  stable  as  the  10  used  in  our  experiment. 
The  QPN  limit  is  shown  under  the  same  conditions  (black  dashed  line),  assuming  50,000  atoms.  The  inset 
illustrates  an  interleaved  interrogation  of  two  atomic  systems,  allowing  continuous  monitoring  of  the  LO 
for  suppression  of  the  Dick  effect.  Dead  times  from  atomic  preparation  or  readout  in  one  system  are 
synchronized  with  dock  interrogation  in  the  second  system.  The  solid  blue  lineindicates  the  suppressed  Dick 
instability  in  the  interleaved  interrogation  scheme  using  Ramsey  spectroscopy  with  an  unimproved  LO. 
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ultran  arrow  band  electronic  clock  transition.  These 
atoms  are  interrogated  simultaneously,  improv 
mg  the  atomic  detection  signal  to  noise  ratio  and, 
thus,  instability,  which  is  limited  fundamentally 
by  quantum  projection  noise  (QPN)  (II),  Like 
other  cycled  atomic  clocks,  the  lattice  clock  suf 
fers  from  a  technical  noise  source  known  as  the 
Dick  effect,  arising  when  an  oscillator's  noise  is 
periodically  sampled  (12,  13 ).  Cavity  stabilized 
lasers  with  k>w  thermal  noise  have  reduced  die 
Dick  effect,  enabling  clock  instability  below 
10” 15  at  short  times  (14),  Recently,  an  uncone 
lated  comparison  of  two  strontium  lattice  clocks 
rewaled  clock  instability  of  3  x  10  16  short 
times*  averaging  to  1  x  1 0  17  in  1 000  s  (5)  or,  in 
another  case,  reaching  the  10  17  lewl  in  20,000  s 
(9).  Here,  by  comparing  two  independent  op 
tical  lattice  clocks  using  ultracold  1  '^Yb,  we 
demonstrate  a  dock  instability  of  1,6  x  10_1S  in 
25,000  s. 

Both  clock  systems,  referred  to  here  as  Yb  1 
and  Yb  2,  independently  cool  and  collect  mYb 
atoms  from  thermal  beams  into  magneto  optical 
traps  [see  Fig.  I  and  (/5)].  Two  stages  of  laser 
cooling,  first  on  the  strong  1Sq  lP\  cycling  trail 
sition  at  399  nm,  followed  by  the  weaker  1Sq  ~P\ 
intereombinatson  transition  at  556  nm,  reduce  the 
atomic  temperature  from  800  K  to  10  pK.  Each 
cold  atom  sample  is  then  loaded  into  an  optical 
lattice  with  -300£r  trap  depth  (Et/kB  1 00  nK; 

recoil  energy,  Etoltzmann's  constant)  foamed 
by  retroreflecting  -6G0  mW  of  laser  power,  fixed 
at  ftie  k(magic'T  wavelength  ^  ~  759  nm  (16)  by 
a  reference  cavity.  At  both  electronic  states 
of  the  clock  transition  are  Stark  shifted  equally 
(/  7,  IE),  For  the  measurements  described  here, 
-5000  atoms  captured  by  each  lattice  are  then 
optically  pumped  to  one  of  the  two  ground  state 
spin  projections  m¥  ±1/2  using  the  ]Sq  3/?1 
transition.  After  this  state  preparation,  applying 
a  140  ms  long  rr  pulse  of  578  nm  light  resonant 
with  the  1  So  3Po  clock  transition  yields  the  spec 
troscopic  line  shape  shown  in  Fig  IC,  with  a 
Fourier  limited  linewidth  of  6  Hz,  Experimental 
clock  cycles  alternately  interrogate  both  mr  spin 
states  canceling  first  order  Zeeman  and  vector 
Stark  shifts.  The  optical  local  oscillator  (LO)  is 
an  ultrastable  laser  servo  locked  to  a  high  finesse 
optical  cavity  (14)  and  is  shared  by  both  Yb  sys 
terns.  Light  is  frequency  shifted  into  resonance 
with  the  clock  transition  of  each  atomic  system  by 
independent  acousto  optic  modulators  (AO Ms), 
Resonance  is  detected  by  mmkoring  file  1S:>  ground 
state  population  (TVg)  and  3/*0  excited  state  pop 
ulation  (Ne).  A  laser  cycles  ground  state  atoms 
on  the  %  1/,1  transition  while  a  photomultiplier 
tube  (PMT)  collects  the  fluorescence,  giving  a 
measure  of  N&.  After  5  to  10  ms  of  cycling,  these 
atoms  are  laser  heated  out  of  the  lattice.  At  this 
point,  Ne  is  optically  pumped  to  the  lowest  lying 
3D]  state,  which  decays  back  to  the  ground  state. 

The  1  jj,  ]P]  transition  is  cycled  again,  now  mea 
suring  Ne,  Combining  these  measurements  yields  a 
normalized  mean  excitation  Ne/(Ne  +  Ay,  Dur¬ 
ing  spectroscopy,  special  attention  was  paid  to 
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eliminating  residual  Stark  shifts  stemming  from 
amplified  spontaneous  emission  of  the  lattice  la 
sers,  eliminating  residual  Doppler  effects  from  me 
chanical  vibrations  of  the  apparatus  correlated  with 
the  experimental  cycle,  and  controlling  the  cold 
collision  shift  due  to  atomic  interactions  within 
the  lattice  (19). 

By  measuring  the  normalized  excitation  while 
modulating  the  clock  laser  frequency  by  ±3  Hz, 
an  error  signal  is  computed  for  each  Yb  system. 
Subsequently,  independent  microprocessors  pro 
vide  a  digital  frequency  correction  f  2(t)  at  time 
t  to  their  respective  AOMs,  thereby  maintaining 
resonance  on  the  line  center.  In  this  way,  though 
derived  from  the  same  LO,  the  individual  laser 
frequencies  for  Yb  1  and  Yb  2  are  decoupled 
and  are  instead  determined  by  their  respective 
atomic  samples  (for  all  but  the  shortest  time 
scales).  The  frequency  correction  signals  f,2(t) 
are  shown  in  Fig.  2A  for  a  5000  s  interval. 
Because  the  experimental  cycles  for  each  clock 
system  are  unsynchronized  and  have  different 
durations,  f\^(t)  signals  are  interpolated  to  a 
common  time  base  and  subtracted  to  compute 
the  frequency  difference  between  Yb  1  and  Yb  2, 
as  shown  in  Fig.  2,  B  to  D.  Measurements  such  as 
these  were  repeated  several  times  for  intervals  of 
-15,000  s,  demonstrating  a  clock  instability  reach 
ing  4  x  1018  at  7500  s.  While  collecting  data  over 
a  90,000  s  interval,  we  observed  the  instability  curve 
in  Fig.  3,  shown  as  the  total  Allan  deviation  for  a 
single  Yb  clock.  Before  data  analysis,  -25%  of 
the  attempted  measurement  time  was  excluded 
due  to  laser  unlocks  and  auxiliary  servo  failures 
(75).  Each  clock  servo  had  an  attack  time  of  a  few 
seconds,  evidenced  by  the  instability  bump  near 
3  s.  At  averaging  times  x  1  to  5  s,  the  instability 
is  comparable  to  previous  measurements  (14)  of 
the  free  running  laser  system,  and  at  long  times, 
the  instability  averages  down  like  white  frequen 
cy  noise  as  -3.2  x  10  16 /y/i  (for  x  in  seconds), 
reaching  1.6  x  10-18  at  25,000  s. 

Also  shown  in  Fig.  3  is  an  estimate  of  the 
combined  instability  contribution  (blue  dashed 
line)  from  the  Dick  effect  and  QPN  (shaded  re 
gion  denotes  the  uncertainty  of  the  estimate). 
These  contributions  must  be  reduced  if  1018  in 
stability  is  desired  at  time  scales  under  100  s. 
Substantial  reductions  of  QPN  are  possible  with 
the  use  of  higher  atom  numbers  and  longer  inter 
rogation  times.  Further  stabilization  of  the  optical 
LO  will  continue  to  reduce  the  Dick  effect,  by 
both  lowering  the  laser  frequency  noise  down 
converted  in  the  Dick  process  and  allowing  in 
creased  spectroscopy  times  for  higher  duty  cycles. 
Improved  LOs  will  use  optical  cavities  exhibiting 
reduced  Brownian  thermal  mechanical  noise  by 
exploiting  cryogenic  operation  (20),  crystalline  op 
tical  coatings  (21),  longer  cavities,  or  other  tech 
niques  (14).  Figure  4  demonstrates  the  advantage 
of  improving  the  present  LO,  with  four  times  less 
laser  frequency  noise  and  four  times  longer  in 
terrogation  time  (corresponding  to  a  short  term 
laser  instability  of  -5  x  1017).  The  red  dotted 
line  gives  the  Dick  instability,  whereas  the  black 


dashed  line  indicates  the  QPN  limit,  assuming  a 
moderate  50,000  atom  number. 

Noting  that  the  calculated  Dick  effect  re 
mains  several  times  higher  than  the  QPN  limit, 
we  consider  an  alternative  idea  first  proposed 
for  microwave  ion  clocks:  interleaved  interroga 
tion  of  two  atomic  systems  (13).  By  monitoring 
the  LO  laser  frequency  at  all  times  with  inter 
leaved  atomic  systems,  the  aliasing  problem  at 
the  heart  of  the  Dick  effect  can  be  highly  sup 
pressed.  The  solid  blue  line  in  Fig.  4  illustrates 
the  potential  of  a  simple  interleaved  clock  inter 
rogation  using  Ramsey  spectroscopy.  Even  with 
the  present  LO,  the  Dick  effect  is  decreased 
well  below  the  QPN  limit  afforded  by  a  much 
improved  LO  (black  dashed  line).  In  this  case, 
spin  squeezing  of  the  atomic  sample  could  re 
duce  the  final  instability  beyond  the  standard 
quantum  limit  set  by  QPN  (22).  The  two  system, 
interleaved  technique  requires  spectroscopy  times 
that  last  one  half  or  more  of  the  total  experimental 
cycle.  By  extending  the  Yb  clock  Ramsey  spec 
troscopy  time  to  >250  ms,  we  achieved  a  50% 
duty  cycle  for  each  system,  demonstrating  the 
feasibility  of  this  technique.  Further  suppression 
of  the  Dick  effect  can  be  achieved  with  the  use 
of  a  more  selective  interleaving  scheme.  Duty 
cycles  >50%  can  also  be  realized  with  the  aid  of 
nondestructive  state  detection  (23). 

Another  important  property  of  a  clock  is  its 
accuracy,  which  results  from  uncertainty  in  sys 
tematic  effects  that  alter  a  standard’s  periodicity 
from  its  natural,  unperturbed  state.  In  2009,  we 
completed  a  systematic  analysis  of  Yb  1  at 
the  3  x  10-16  uncertainty  level  (16).  Since  then, 
we  reduced  the  dominant  uncertainty  due  to  the 
blackbody  Stark  effect  by  one  order  of  mag 
nitude  (24).  With  its  recent  construction,  Yb  2 
has  not  yet  been  systematically  evaluated.  The 
fact  that  the  instability  reaches  the  1018  level 
indicates  that  key  systematic  effects  (e.g.,  the 
blackbody  Stark  effect,  atomic  collisions,  lattice 
light  shifts)  on  each  system  are  well  controlled 
over  the  relevant  time  scales.  The  mean  frequency 
difference  in  Fig.  2  was  (f2(t)  f(t))  30  mHz, 

which  is  within  the  Yb  1  uncertainty  at  1016  Sys 
tematic  effects  on  each  system  can  now  be  effi 
ciently  characterized  beyond  the  1017  level.  With 
continued  progress,  we  envision  1018  instabil 
ity  in  only  100  s  and  long  term  instability  well 
below  10-18. 

Clock  measurement  at  the  1018  level  can  be 
used  to  resolve  spatial  and  temporal  fluctuations 
equivalent  to  1  cm  of  elevation  in  Earth’s  gravi 
tational  field  (25  28),  offering  a  new  tool  for  ge 
odesy,  hydrology,  geology,  and  climate  change 
studies.  Space  based  implementations  can  probe 
alternative  gravitational  theories,  e.g.,  by  mea 
suring  red  shift  deviations  from  general  relativity 
with  a  precision  that  is  three  orders  of  magni 
tude  higher  than  the  present  level  (28).  Though 
present  day  temporal  and  spatial  variation  of 
fundamental  constants  is  known  to  be  small 
(6,  29,  30),  1018  level  clock  measurements  can 
offer  tighter  constraints.  Finally,  timekeeping 


improvements  can  benefit  navigation  systems, 
telescope  array  synchronization  (e.g.,  very  long 
baseline  interferometry),  secure  communication, 
and  interferometry  and  can  possibly  lead  to  a  re 
definition  of  the  SI  second  (9). 
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